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Abstract 

We consider a microscopic model for the failure of soft adhesives in tension based on ideas 
of bond rupture under dynamic loading. Focusing on adhesive failure under loading at constant 
velocity, we demonstrate that bi-modal curves of stress against strain may occur due to effects 
of finite polymer chain or bond length and characterise the loading conditions under which such 
bi-modal behaviour is observed. The results of this analysis are in qualitative agreement with 
experiments performed on unconfined adhesives in which failure does not occur by cavitation. 



Introduction 



The failure of natural and artificial adhesives is a subject of great importance involving the con- 
fluence of physical chemistry as well as statistical and continuum mechanics, but is as yet poorly 
understood. This is because of the multiple hierarchy of spatial and temporal scales involved, from 
the failure of individual molecular bonds to that of the macroscopic sample, that occur over tempo- 
ral scales that can also span many orders of magnitude. From a practical viewpoint, a particularly 
simple, yet useful, method of probing the failure of an adhesive is to stretch it at constant speed 
and plot the nominal stress against strain (referred to as tack curves) to deduce some insight into 
the physical processes involved in failure. In fig. ^ we show some typical tack curves observed 
during the loading of an adhesive in an unconfined geometry [3] and in a confined geometry [3J- 
We see that typical tack curves are approximately linear at small strains (since adhesives respond 
approximately linearly to small deformations), but at larger strains the nominal stress peaks (as the 
adhesive becomes damaged and is less able to resist deformation) and falls to zero as the adhesive 
fails. As the parameters associated with the tack test, such as the rate of loading, the degree of 
confinement (i.e. the aspect ratio of the adhesive) or the type of adhesive are varied, a qualitative 
change in the tack curves shown in fig Q may be observed in the appearance of a plateau or even 
a second peak. Since the area under a tack curve is a measure of the toughness of the adhesive, 
understanding the processes involved in determining the qualitative transitions in the stress-strain 
curve such as the maximum stress, and the presence or absence of a second peak or plateau are 
clearly efficacious in designing systems with a prescribed toughness. 



1 




Figure 1: Some typical tack curves showing both bi-modal and uni- modal behaviours. In both 
figures, the change from uni-modal to bi-modal is induced by changing the velocity of loading, (a) 
Loading force, Ptensiie, as a function of imposed displacement, 5 tens ii e , in the case of an unconfined 
adhesive (aspect ratio of order 1), taken from Crosby and Shull [3]. (b) Nominal stress, a, as a 
function of strain, e, for a confined adhesive (a very thin layer of adhesive), taken from Chiche 13 . 



Some previous work has focused on explaining the appearance of a plateau in the tack curves 
(fig- IHb)) as being due to the appearance of cavitation bubbles within the bulk of a confined 
adhesive layer (see, for example, 0). However this picture does not account for the subsequent 
rise in force that is observed even in the absence of cavitation. An alternative view focuses on 
the implications of non-linear elasticity on the shape of tack curves. For example, Sizaire and 
Legat |15j used numerical computations based on the FENE-CR constitutive relation to model the 
response of a viscous suspension of finite length polymer chains. More recently, Roos and Creton 
j!2j used experimental data to fit a non-linear constitutive equation based on the 'slip-tube' model 
of polymer elasticity to data obtained via tack tests. While some of the qualitative features of 
the tack curves shown in fig. n are reproduced by these models, neither accounts for the failure of 
the material via a reasonable description in terms of the rupture of bonds or equivalent molecular 
picture. 

Here we attempt to remedy this omission by considering these molecular effects within a simple 
microscopic model of the failure process. In particular, we show how changes in the nature of the 
adhesive and the rate of loading may lead to bi-modal tack curves whose physical origin is the 
finite length of either adhesive bonds or the polymers within the adhesive. We also characterize 
the loading conditions under which these bi-modal curves might arise and show that these are 
consistent with experimental conditions. Although our simple model neglects many effects that are 
undoubtedly of importance, it sheds some light on the underlying physics and suggests experiments 
that might further elucidate both the understanding and tailoring of adhesive toughness. 
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Figure 2: Schematic of the model used in this discussion (adapted from Seifert \1A ]). Two rigid 
plates are pulled apart at constant velocity v so that the loading of the springs connecting them is 
under displacement control. The bonds dissociate at a rate dependent on the load applied to them 
but are not allowed to reform. 

2 A Microscopic Model 



The idea of bond rupture enhanced by the presence of an applied load was first used in the polymer 
literature to explain the frictional properties of elastomers by Schallamach JHj following the classical 
work of Eyring [3] and Kramers [5] on reaction rate theory. Similar ideas have also been used in 
more recent work on adhesion [21 |Hj to investigate the rate dependence of adhesive properties such 
as the work of adhesion, and also to understand the molecular and cellular aspects of biological 
adhesion , ^3] . A reduced model based on such ideas does not attempt to include the intricacies 
of the experimental geometry but allows us to focus on the physical processes that might lead to 
bi-modal tack curves. 

Our starting point is to consider two rigid plates connected by N\,{t) non-linear springs of natural 
length L (see fig. |2J). Chaudhury [2] and Ghatak et al. |S] have considered a similar model in 
which the springs represent the adhesive bonds that join the polymer adhesive to the plates. This 
model of adhesive failure (failure at the interface between the adhesive and apparatus) could, 
however, be applied more generally to cohesive failure (failure within the bulk of the adhesive) 
when the cohesive failure occurs via the rupture of polymer chains within the bulk rather than the 
appearance of cavities. We shall not, however, distinguish between these two modes of failure, since 
both can be represented by the schematic picture of fig. [21 

Once the plates begin to be separated, the springs extend and so are stressed, causing them to 
dissociate at a rate that is dependent on the load to which they are subjected. The dissociation 
rate is enhanced due to the effective lowering of the bond-potential barrier brought about by the 
application of a force. Assuming a simple Arrhenius type dependence, which is often referred to as 
the Bell model l s , we write the dissociation rate as: 

where 7 is a lengthscale associated with the bond potential, ksT sets the energy scale, F(t) is the 
force with which each bond is loaded (assuming that they are all equally stressed) and k°_ is the 
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dissociation rate in the absence of a force, which itself can be estimated in terms of the Eyring or 
Kramers jH] theory of chemical reactions. Given a particular constitutive relation for the non-linear 
springs it is a simple matter to write the corresponding evolution equation for the number of intact 
springs, Nj,, remaining at time t. However, before doing so, we need a model for the response 
of the individual springs or polymer chains at the interface or in the bulk of the adhesive. We 
choose a form that is simple for the purposes of analysis, but incorporates the two features that we 
desire here: linear behaviour at small displacements and a simple divergence in the force required 
to stretch the polymer as its length approaches the natural length, L. One possibility is the choice 
of a constitutive relation similar to the freely jointed chain model for stretching of a polymer chain 
given by the inverse Langevin function (see, for example, ^Ll\), although this has no closed form. 
Therefore, for simplicity, we assume that the force required to stretch each spring by a distance x 
is F = kBTx/b(L — x), where L is some natural length, such as the contour length of the polymer 
chain between crosslinks, 6 is a molecular length associated with the bond/polymer chain such as 
the monomer or blob size jllj and k^T is the thermal energy scale. The typical stiffness here 
has been chosen to reproduce that of a Gaussian chain at small strains. In fact, the above 
force-displacement relation is similar to a simple Pade-like interpolation of the linear behavior for 
small displacements and the divergence implicit in the inverse Langevin function |lflj . 

Assuming that the loading of bonds/polymer chains takes place at the same constant velocity as 
the plates are separated, v, the constitutive law described above allows the force on each remaining 
bond to be written as a function of time, t: 



Here we have assumed that the rate of reforming of bonds is relatively small, which is certainly 
true if the pulling force is moderate or large, the case of most interest here. There are three 
natural timescales in this problem. The first of the time-scales is the mean lifetime of an unforced 
spring, namely l/k_. The second characterises the rate at which energy is stored in a spring 
compared to thermal energy bL/jv and provides us with a useful non-dimensionalisation of time: 
r = ■jvt/bL. The third is simply the time taken to stretch a constituent spring to its natural length 
L, given by L/v. From the ratios of these time-scales, we find two dimensionless parameters that 
govern the evolution of the system. The first of these is e = jv/k^bL, the ratio of the intrinsic 
off rate to the loading rate, which is readily varied by varying the rate at which loading takes 
place, v. The second dimensionless parameter is the time at which the natural length is reached, 
T m = 7/6, and is the ratio of the width of the bond-potential barrier and the monomer length. 
As this is dependent solely on the chemical properties of the adhesive, it is more difficult to vary 
experimentally, though this may, in principle, be achieved by altering the typical monomer bond 
length, 6, and the strength of the monomer-monomer interaction; large values of r m correspond 
to very soft (long range) interactions between small monomers, such as Coulomb-type interactions 
between nano-particles, while small values of r m correspond to short range interactions between 
large monomers, such as van der Waals interactions between macromolecular globules. Making use 
of JJ) and (J2J) the governing equation for the number of intact bonds, i.e. dN^/dt = — A;_(i)A^, may 
be written in dimensionless form as 



F(t) 



ksT vt 
bL 1-f 





(3) 
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Figure 3: Typical dimensionless force-displacement graphs for the model shown in Figure 2. As 
the rate of loading is increased (in the sense indicated by the arrow), the tack curves change 
from unimodal to bimodal and then back to unimodal. Here r m = 0.05 and the rate of loading 
is proportional to the parameter e, which takes values 0.006 (unimodal with a slow force decay 
preceding failure), 0.01 (bimodal) and 0.013 (unimodal with a very sharp force drop preceding 
failure) . 



The change of variable u = t/(1 — r/r m ) allows us to write the solution of © in the form 

Wt) = AWO) exp U f "> 7 ( 7> *) , (4) 

y e Jo [1+W/TmY J 

which can then be expressed in terms of the exponential integral Ei(z), although this representation 
of the solution is of little interest in what follows. Here -/V&(0) = Nq, the initial number of intact 
bonds; however, we take Nq = 1 so that Nj,(t) represents the fraction of the original bonds remaining 
at time r. Then we may write the total dimensionless force on the plates J~{t) required to maintain 
the constant velocity displacement (scaled by j/ksT) as the product of the force per bond and the 
number of intact bonds so that 

Hr) = : T , N h {r). (5) 

1 - T/T m 



As a function of time (or equivalently of displacement) the force (J2J) is not monotonic but instead 
rises slowly and then decays down to zero. An explanation for an intermediate maximum follows 
from the simple fact that initially the force per bond increases with the displacement so that the 
total force also increases. However, as individual bonds start to break in response to the applied 
load, the force eventually starts to decrease. What is more interesting is that for certain values 
of the parameters e and r m , it is possible for there to be two maxima in this curve. Examples in 
which there are one or two maxima in the force are given in fig. |21 



Having observed that it is possible to obtain bi-modal tack curves that bear some qualitative 
resemblance to those shown in fig.n( a )> we now focus on understanding such curves in more detail. 
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Figure 4: Phase diagram showing the values of the parameters r m , the scaled bond-potential width 
and e, the scaled loading rate for which the scaled applied force, J-(t), is unimodal or bimodal. As 
expected, for very large or very small loading rates, the tack curve is unimodal, but for a narrow 
range of loading rates, bimodal tack curves can exist. 



3 Phase space for bi-modal tack curves 



Next, we quantify the regions of (r m , e) space in which bi-modal tack curves occur and then explain 
the physical origin of this behaviour. We first obtain an expression for the (non-dimensional) times, 
t*, at which J~(t) has a turning point. Differentiating the expression for J-(t) in © and using (j3J) 
to eliminate dN^/dr we see that the r* satisfy 

e = r*(l-r7r m )expf T \ ) = G(r m ,r*). (6) 

VI - T*/T m J 

Although this equation does not admit an analytic solution, the general properties of the function 
G(T m ,r) are enough for our purposes. Clearly G(r m ,0) = and as r — > r m , the value of G(r m ,r) 
diverges exponentially. However, in the interval (0, r m ), G may have, depending on the value of r m , 
either zero or two turning points. Assuming that there are two turning points, they occur when 
t = t± where 

3T m T m i T m \J 't^j 6T m + 1 . 

r± = ^ • (7) 

For the values of t± to be real and distinct we require that — 6r m + 1 > 0, i.e. r m < 3 — 2y/2 
or r m > 3 + 2y/2. For values of r m satisfying the latter condition, it is easy to see that r + < 
which is not physically meaningful so that we consider only the first possibility. In this region 
< t± < T m so that we really do see both of the turning points in the function G(r m ,r) in the 
interval [0,r m ). From the geometrical constraints on G, we see immediately that r_ corresponds 
to a local maximum and r + a local minimum in G(r m , r). Thus we see that the applied force, 
is bimodal only if e E (e_, e+) where e± = G(r m , tv). 

This result is most easily summarized by means of a phase diagram in (r m , e) phase space showing 
where we expect the two types of force curves (unimodal and bimodal) to occur. The borders of 
these regions can be calculated analytically from the preceding analysis and are shown in fig. 21 
which shows that there is only a very small range of the parameters r m and e for which the bimodal 
behaviour can be observed. 
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Comparison with experiments requires estimates for the typical values of the parameters e and r m . 
Assuming that the failure occurs as a result of the rupture of individual polymer chains within the 
bulk rather than adhesive bond rupture at the interface we use the result for a Gaussian chain 11 
which has the constitutive response given by (j2J with L = Nb where TV" ~ 10 8 is the number of 
monomers in an individual chain and b ~ 10 _9 m is the typical monomer (or Kuhn) length. Along 
with 7 ~ 10" 10 m, this gives r m ~ 0.1. It is difficult to estimate a typical value of e because of the 
ambiguity in k_. However, it is enough to note that whatever the value of k_, the typical value 
of T m is consistent with our model's prediction that the appearance of bimodal tack curves can be 
controlled by varying the velocity of pulling, as is observed experimentally. 

To understand the behaviour exhibited in fig. 0] we first consider the existence of an upper bound 
on the value of r m (= "f/b) for which bimodal force displacement curves may occur. For larger 
values of r m , the nonlinearity in the force law becomes important only once most of the bonds 
have already been broken and so no second peak is observed. For a given value of r m (sufficiently 
small to avoid the previous pitfall), very fast loading (large e = jv/k°_bL) will cause most of the 
bond breakage to occur only as the strain-stiffening starts to become important and, again, only 
one force peak is observed. We thus expect that bi-modal tack curves can only occur below some 
critical loading rate e+. Notice that the bi-modal curves presented in the experimental data fig. 
n are not observed for large loading rate: just as predicted by our analysis and simple physical 
mechanism. On the other hand, if we load the springs too slowly then very few of them will be 
broken by being stressed and instead break over their natural time-scale l/k^_ so that few remain 
by the time they are extended to the point at which strain hardening is appreciable. This leads us 
to expect that for a given value of t there should be a minimum loading rate e_ below which the 
force curve is again unimodal. 

As we have just seen, bimodal curves occur when the individual bond-springs are extended towards 
their natural length (when r = r m ) and stiffen as a result. We observe that the decay of the applied 
force after this peak is much more rapid than the decay after the first peak. Informal comparison 
with the experimental results of Crosby and Shull [3] presented in fig. ^ (a) confirms that this is 
also observed experimentally. However, the results of Chiche 3 show that in a confined geometry 
where cavitation occurs, such a rapid decay is not experienced after the second peak. A natural 
question that arises is whether the rapid decay is a result of a particular choice of parameters or is 
an intrinsic characteristic of such a model. In the next section, we show that our model predicts 
that this decay in the force is necessarily very rapid. 



4 Late time asymptotics 

To study the behaviour of the applied force after the second peak, we study the behaviour of Tir) 
as r — > r m , i.e. for r such that r m -T<l. In terms of the parameter u = r m /(l — r/r m ) this is 
equivalent to the limit u> 1. With this substitution, we can integrate © so that, up to a constant 
of integration, the number of bonds Nb satisfies 




(8) 
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which can then be integrated by parts to give 

e « z V u / 

Keeping only the leading order terms, recalling that u> 1, yields an expression for the force, which 
to leading order is 

Hr) ~ r^ exp (-^=^(1 -r/r m ) 2 exp ( Tm )) . (10) 
1-T/T m V e \(1 — r /T m ) J J 

It is immediately clear from this expression that when t ~ Tm? 

the pulling force decays super- 

exponentially and that the initial conditions are irrelevant since their influence is exponentially 
small. In fact, the decay given in H10|) is impossible to observe in the numerical solutions of (J3J) 
since it is important only very close to the singularity, where the errors in the numerical integration 
scheme overwhelm the effect that we wish to observe. 



This result shows that the behaviour of the force-displacement curve as r — ► r m is qualitatively 
similar to that observed in some experiments (e.g. [1]), but not others (e.g. [3]), which differ from 
each other mainly in the nature of the experimental geometry highlighting the role of confinement. 
In an unconfined adhesive, the failure is initiated by a fingering instability that occurs around the 
edge of the adhesive [3]; a confined adhesive will fail via cavitation within the bulk of the adhesive 
In the former case, the constitutive response of the elastomer that remains is unchanged but 
in the latter case, the properties of the adhesive change because of its ability to accommodate 
the imposed deformation via bubble growth. Thus our microscopic model applies only to adhesive 
loading in the unconfined case. In this scenario, damage is accumulated by the adhesive as a whole 
even though the individual constituents remain largely undamaged. This notion of damage to the 
adhesive as a whole can be used to recast the microscopic model of §2 in terms of a continuum 
theory, to which we turn next. 



5 A Continuum Model 



To adapt the notion of adhesive failure in terms of the microscopic spring model of ^Uto the case of 
a continuous adhesive of constant volume we introduce a damage parameter, 1 — a, which increases 
from to 1 as damage is accumulated by the adhesive. Here a is defined as the ratio of the actual 
stress a required to impose a given displacement and the stress experienced by an undamaged 
adhesive subject to the same displacement, a u i.e. a = aa v 



By analogy with the parallel spring model of we postulate an Arrhenius type evolution for a as 
a function of the imposed stretch, A, namely 

da a ( a u V\ 

whose basis is the correspondence between a and the number of unbroken springs iVj, in the micro- 
scopic model. Here, uj = \/k®__ is the ratio of the stretch rate to the unforced dissociation frequency, 
the continuum analog of e, and V ~ 7^ is some characteristic volume (where 7 is again a typical 
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Figure 5: The scaled nominal stress N as a function of the applied stretch A obtained from the 
continuum model (J11I13|) . (a) Increasing the loading rate, parametrized by to, changes the behaviour 
from unimodal to bimodal and back again. Here J m = 1, v = 0.05 and uj takes the values 
0.1,0.2,0.23 and 0.26. (b) Numerically computed phase diagram showing the values of the scaled 
loading rate u> and the value of the first invariant at failure, J m , for which bimodal and unimodal 
behaviours are observed. Here the scaled bond-potential barrier v ~ j/l e = 0.05. 



length associated with the bond potential and l e is the typical entanglement length or distance be- 
tween crosslinks of the polymer network, so that stress is distributed over an area 1%). To complete 
the formulation of the continuum model, we need a constitutive equation for the stress in terms of 
the stretch or equivalently a free energy Uf(X) such that 

- ^ < 12 > 

To account for the finite length effects of the constituent chains in the adhesive, we use the form 
of the free energy, Uj, suggested by Gent (1996) 

U f = ^J m log(l-Pi, (13) 



' rn 



where V ~ l\ is the typical microscopic volume characteristic of the polymer network, J(X X , X y , X z ) = 
A 2 , + Xy + X 2 Z — 3 (subject to the constraint X x X y X z = 1 imposed by incompressibility) is the first 
invariant of the deformation gradient and J m = J(A m ), for some A m , is the maximum value of this 
invariant corresponding to the maximal possible stretch of the adhesive, X m . We note that when 
J ^ Jm, the above free energy reduces to Uf ~ ^yf-J, the classical expression for a neo-Hookean 
solid. For an incompressible elastomer ic adhesive loaded axially, as we will assume here, X x = X y 
so that J = J(A) = A 2 + 2A" 1 - 3. 

In this formalism, it is a simple matter to solve (|11II13|) subject to the initial condition that a(A = 
1) = 1 via numerical methods. Such a solution can then be used to infer the nominal stress, J\f, 
given by 

Experimentally, this is the quantity most commonly measured and thus allows us to make some 
simple qualitative comparisons. In figEl (a) we show representative solutions of ()11II14|) for pa- 
rameter values J m = 1 and v = V/V ~ j/l e = 0.05 as the loading rate is increased. We note 
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that the macroscopic analog of the bond potential barrier v is characterized as a fraction of the 
entanglement length rather than the monomer length as in §2. With this caveat, we see that our 
results are qualitatively very similar to those obtained using the microscopic model in §2; in par- 
ticular, there is a range of loading rates for which the nominal stress is bimodal, and for very large 
loading rates the nominal stress decays very rapidly after the second peak. It is also possible to 
characterize the region of (J m ,u,u) space that leads to bi-modal tack curves, although this must 
be done numerically here. In fig. |S] (b) we characterize this region and see results that are, again, 
qualitatively similar to those in §2. 



6 Conclusions 



The existence of a second peak in tack curves obtained in adhesive testing and failure is a common 
feature of many recent experiments. When combined with the observed reduction of the nominal 
cross-sectional area of the adhesive at the point where the second peak occurs this indicates that 
the modulus of the adhesive as a whole must have increased by this point. In this paper, we have 
investigated the possibility that this might be produced by the effects of finite bond/chain length 
and shown that such finite length effects, in the absence of cavitation, produce tack curves in 
qualitative agreement with those obtained experimentally with unconfined adhesives (e.g. Crosby 
and Shull 4 ). Our model allows us to also determine the region in phase space where this bi-modal 
behavior may be seen, and further explains the rapid failure of the adhesive following the second 
peak in the force. We predict that these bi-modal features should disappear at high pulling velocity, 
as is observed in such experiments. 

One deficiency in the model considered here is its exclusion of the flow-induced viscoelastic nature of 
the adhesive itself; for simplicity we have considered purely elastic effects, with the rate dependence 
coming in via a simple kinetic mechanism. This might simply be remedied by incorporating the 
viscous and elastic response of the polymer chain with a bond joining this chain to the substrate. 
Another improvement to the model presented here would be to take into account the spatially 
varying nature of the problem. There are several levels of sophistication at which such an approach 
could be introduced. Perhaps the simplest modification would be to consider an axially symmetric 
problem but to account for the variations in chain length within the adhesive that occur because of 
the necking phenomenon that is a result of incompressibility of the adhesive and the pinning of the 
contact region to the plate. Ultimately, with the full spatial dependence accounted for, it would be 
possible to investigate the interplay between the onset of the fingering instability and finite chain 
length effects. Additionally an important aspect of any continuum theory of failure that we have 
not accounted for is the role of stress concentrations in the vicinity of the crack that demarcates 
the failure region. This effect is important in brittle materials, but its relative role in elastomers is 
debatable since the locally high stresses lead to crack blunting and even slow flow in soft materials. 
In the absence of experimental studies of the relative importance of crack blunting, we defer this 
important theoretical question to the future. 

Our work has focused on one possible mechanism for the appearance of bi-modal tack curves. Other 
effects that might enhance this could arise from the realignment of polymer domains within the 
adhesive: at intermediate displacements, it is energetically cheaper to accommodate the imposed 
deformation by rotation of polymer domains than to stretch them. This semi-soft response (see, for 
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example Warner and Terentjev ^H] f° r a thorough discussion of semi-softness) alters the constitu- 
tive response of the material and would allow for a relaxation of nominal stress via domain rotation. 
Once the maximal rotation has occurred, the domains must once again deform allowing them to 
be stretched as springs before they eventually break thence giving rise to a second peak. Indeed, 
this suggests that experiments that probe the structural rearrangements of the adhesive (e.g. us- 
ing birefringence measurements) simultaneously during mechanical testing might help unravel the 
relative contributions of domain rotation and simple stretching mechanisms. 
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